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t*"* ■ In this paper, bound states energies and corresponding wave functions of H-shaped quan- 

tum wires are calculated numerically in the presence of external magnetic and electric fields 
and within the Landau gauge. With a suitable definition of external confinement potential, 
we present a numerical algorithm to calculate the profile of probability distribution of charge 
carriers. Our analysis shows that in the presence of external electric and magnetic fields, bound 
state properties of carriers are sensitive functions of an asymmetric parameter a = ^ which 
measures the relative width of the well in two directions. We also study many body effect of 
bandgap renormalization in this quasi one dimensional system within dynamical random phase 
approximation . 
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1 Introduction 



A highly dense electron-hole plasma can be generated in a wide variety of semiconductors by optical 
pumping. The band structure and the optical properties of highly excited semiconductors generally 
differ from those calculated for non-interacting electron-hole pairs due to many-body exchange- 
correlation effects arising from the electron-hole plasma[l,2]. In recent years, quasi-one dimensional 
semiconductor quantum wires (QW) have been fabricated in a variety of geometric shapes with 
atomic scale definition, and QW optical properties have been studied for their potential device ap- 
plications such as semiconductor lasers [3 ,4]. Recently different geometries of quantum wires, such as 
rectangular, V-shaped, L-shaped and T-shaped quantum wires have been fabricated and some of their 
electronic properties have been studied[4]. In addition, various experimental techniques for fabrica- 
tion and growth of these structures have been developed [3-5]. Square quantum well wires have been 
studied by Hu and Das Sarma[6]. They have calculated the value of the band gap re-normalization( 
renormalization of fundamental gap structure due to many body exchange-correlation effects) for this 
case within GW approximation. Excitonic effects in quantum wires have been studied by Goldoni et 
al [7] . They have studied the effects of Coulomb interaction on the linear and nonlinear optical prop- 
erties of both V-shaped and T-shaped semiconductor quantum wires. Wang and Das Sarma have 
proposed an elegant framework for numerical studies of carrier induced many-body effects on the ex- 
citonic optical properties of highly photoexcited one-dimensional quantum wire systems [8, 9]. Hwang 
and Das Sarma have investigated the dynamical self-energy corrections of electron-hole plasma due 
to electron-electron and electron-phonon interactions at the band edges of a quasi-one dimensional 
photoexcited electron-hole plasma within GW approximation[10,ll]. Rinaldi and Cingolani have 
studied the optical properties of quasi one dimensional quantum structures specially the case of 
V-shaped quantum wires[12]. Bener and Haug have considered plasma-density dependence of the 
optical spectra for quasi-one-dimensional quantum well wires[13]. Tanatar has studied the band gap 
re-normalization in quasi-one dimensional systems in a simple plasmon-pole (quasi-static) approx- 
imation^, 15]. Giiven et al have studied the band-gap renormalization in quantum wire system 
within dynamical correlations and multi-subband effects [16]. Luttinger liquid behavior of a semi- 
conductor quantum wire has been studied by Bellucci and Onorato[17]. They have also studied the 
effects of magnetic field on low dimensional electron systems focusing on Luttinger liquid behavior 
in a quantum wire[18]. In addition, they have studied also the ballistic electron transport in a quan- 
tum wire under the action of an electric field [19] and quenching of the spin Hall effect in ballistic 
nanoj unctions [20]. As an application, recently they have studied the transport through a double 
barrier in large radius Carbon Nanotube with transverse magnetic field [21]. Many particle aspects of 
a semiconductor quantum wire within an improved random phase approximation has been studied by 
Ashraf and Sharma[22]. They have considered structure factor, pair distribution function, screened 
impurity potential and density of screening charge and exchange and screened exchange energies 
within an improved random phase approximation. 
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On the other hand, T-shaped and L-shaped quantum wires have been studied by some authors 
recently. For example, Sedlmaier and his coworkers have studied the band gap re-normalization 
of modulation doped T-shaped quantum wires. They have presented a self-consistent electronic 
structure calculation for this device [23]. Using a density functional theory, Stopa has calculated 
the electronic structure of a modulation doped and gated T-shaped quantum wire [24]. Szymanska 
et al have studied the excitons in T-shaped quantum wires [25]. They have calculated energies and 
oscillator strength for radiative recombination and two particle wave functions for ground state 
exciton in a T-shaped quantum wire. Lin, Chen and Chuu have found the dependence of the bound 
states of L-shaped and T-shaped quantum wires to some asymmetric parameter in an inhomogeneous 
magnetic field[26]. They have proposed a simple model to explain the behavior of the magnetic 
dependence of the bound state energy both in week and strong field regions. Recently Nozari and 
Madadi have studied numerically the bound states properties and the band gap renormalization 
of V-shaped and T-Shaped quantum wires within dynamical random phase approximation [2 7, 2 8]. 
Ballistic transport through coupled T-shaped quantum wires has been studied by Lin et al [29]. 

As another possible geometry of low dimensional systems, H-shaped electronic systems have been 
considered and their many-body electronic properties have been studied by some authors. Shin and 
coworkers have studied quantum transport in an H-shaped quantum wire and a ring structure [30]. 
They have studied numerically the transport properties of an H-shaped quantum wire structure by 
using the mode matching technique. They have reported the existence of anomalous Hall resistance 
plateaus in this structure with relatively low magnetic fields as precursors of integer Hall plateaus. 
Henkiewicz et al have studied the manifestation of the spin Hall effect in a two-dimensional elec- 
tronic system with Rashba spin-orbit coupling via dc-transport measurements in a realistic H-shaped 
structure [31]. Designing of H-shaped micromechanical structure has been studied by Arhaug and 
Soeraasen[32]. Our investigation shows that there is no other concrete study of these structure in 
existing literature. Obviously an analytical-numerical study of this special structure is important to 
fill the existing gap. Specially, our investigation of literature shows that bound states of H-shaped 
quantum wires and many-body effects such as bandgap renormalization of these low dimensional 
systems have not been studied yet. So, in this paper we consider the geometry of H-shaped quantum 
wires and by a suitable analytical definition of quasi-one dimensional H-shaped confinement poten- 
tial, we propose a numerical scheme for calculating the bound states energies and wave functions 
of charge carriers in the presence of external electric and magnetic fields and within the Landau 
gauge. We obtain the profile of charge carriers distribution(probability distribution) in the presence 
of electric and magnetic fields. As an important nonlinear optical effect, the many-body exchange- 
correlation induced band gap narrowing in this type of quantum wire will be studied within leading 
order dynamical random phase approximation. 

The paper is organized as follows: section 2 is devoted to formulation of the problem and definition 
of confinement potential analytically. In section 3 we give a short but complete review of bandgap 
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renormalization in a general quasi-one dimensional semiconductor. Some numerical details are given 
at the end of this section. Section 4 provides numerical results of our study and their interpretation, 
while the numerical scheme of our calculations based on finite difference algorithm is presented in 
the Appendix. The results of each step are shown by figures. Finally, summary and conclusions are 
given in section 5. 



2 The Setup 

The geometry of a typical H-shaped quantum wire is shown in figure 1. The typical values of W x 
and W y are of the order of 5 — 50 Nanometer. We study bound states and many-body effects in this 
quasi one dimensional system in the presence of external electric and magnetic fields. The electric 
field is assumed to be directed along the rr-axis and its typical value is of the order of a few V/ cm. 
The presence of magnetic field is in such a way that the Landau gauge, A = (0, Bx, 0) is satisfied. 
The Hamiltonian of an electron in this configuration can be written as follows 



H =7^ 
2m 



Px + (Py ~ eBx) 2 — eEx (1) 



where e and m are electron charge and mass respectively. The Schrodinger equation of this electron 
with wave function ^(x, y) can be written as follows 

( d 2 d 2 2ieB d\ T/ s ( eB\ 2 , T/ . 2meE T . . 2m x T . . 



dx 2 dy 2 h dy I ' \ h J ' h ti 



where A stands for eigenvalues of energy, E represents the value of electric field and B is the value 
of the magnetic field. To solve this eigenvalue problem we need the boundary conditions which are 
given by the external confinement potential. The geometry of H-shaped quantum wire as shown in 
figure 1, suggests the following definition of confinement potential 



V(x,y) 



-oo < x < +oo, -\W y <y< -W y 

if { -oo < x < +oo, W y < y < \W y 

< x < }ji _ Wy <y<W y 

oo elsewhere 



With this definition of confinement potential (which provides the required boundary conditions), 
we solve equation (2) numerically to find eigenvalues and eigenf unctions of bound states of electron 
in the presence of electric and magnetic fields. Our numerical strategy based on finite difference 
algorithm is presented in the Appendix. The results of these calculations for different external field 
configurations and variety of possibilities are shown in Figures and will be interpreted in section 
4. The bound states wave functions obtained in this section will be used to study the bandgap 
renormalization of H-shaped quantum wire in the next section. 
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3 Many Body Effects in H-Shaped Quantum Wires 



There are several nonlinear optical many body effects in quasi-one dimensional semiconductor systems 
originating from exchange-correlation effects in a dense excited plasma. One of the most important 
many-body effect in high density electron-hole plasma is a density-dependent re-normalization of the 
fundamental band gap of the semiconductor, which causes an increasing absorption in the spectral 
region below the lowest exciton resonance. The exchange-correlation correction of the fundamental 
band gap due to the presence of free carriers (electrons in the conduction band and holes in the valence 
band) in the system is referred to as the band gap re- normalization effect. Optical nonlinearities, 
which are strongly influenced by screened Coulomb interaction in the electron-hole plasma, are 
typically associated with the band gap re-normalization phenomenon. In which follows we use the 
two band(one conduction band and one valance band) model to study the one dimensional electron- 
hole system. We neglect the effects of higher subbands and degeneracies in valance bands. We 
assume that electrons and holes densities are constant in time. In this situation, Hamiltonian of the 
system can be written as [8, 9, 33] 



k I 



+ VcAtiA-qci+qCk'Ck + V cM {q)dl_ q d\, +q d kl d k + V c , eh {q)cl_ q c k d\, +q d k , 



k,k',q 



(4) 



In this equation which contains all information about this one dimensional system, c k and c\ are 
annihilation and creation operators for conduction electrons respectively. Also d k and d\ are anni- 
hilation and creation operators for valance band holes. is the direct band gap between the top 
of the valance band and the bottom of the conduction band. V c>i j show the possible three Coulomb 
interactions between electrons and holes. Note that the two first interactions lead to electron-hole 
quasi particle self-energies while the third one leads to the production of excitonic bound states. Note 
also that this Hamiltonian consists of spin effects, although spin index is not included explicitly. 

The Coulomb interaction matrix element in one dimensional quantum wire is given by the fol- 
lowing relation [8] 

+ °° J J / f + °° J J / f + °° J e-^i^x,^)! 2 !^',?/)! 2 



r+oo ^ r+oo r+oo 

V c ,ij(q) = — dxdx' I dydy' / dz 

£n J —oo J — oo J — oo 



e J-oo J-oo J-oo ^J z 2 _|_ (jj _ y/y _|_ ( x _ x iy 

= _/ dxdx > dydy'^y^^x'.y'^KJq^x - x'f + (y - y')A, (5) 

€q J —oo J -co v ' 

where <f>i(x,y) is the quantum wire confinement wave function for the lowest eigenstate of electrons 
or holes. The exact form of these eignfunctions depends on the geometry and details of confinement 
potential. Kq(x) is the zeroth-order modified Bessel function of second kind. In the setup of our one 
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dimensional quantum system we have assumed that carriers are free to move in z direction but x 
and y are directions of confinement. 

Band gap renormalization in quasi-static approximation is given by [14, 15] 

Ai(fc) = £ [ - V s (k - fc')n i (e iifc ) + \{V s {k') - V c {k'))] (6) 
k' 1 

where 

V s (k) ee V s {k,u = 0) = a (7) 

is the statically screened Coulomb interaction and n^e^k) is the fermion momentum distribution 
function 

miet, fc ) = (e««.*-«> + 1)- 1 (8) 

where e ei fc = + an d = are the bare energies of electron and hole in their respective 
bands and /!« is the chemical potential. e(k,u) as dynamical dielectric function is defined as follows 

-co 2 - [(k 2 /2m t ) - kv F £ 



777 ■ 

e(k, u) = l- V c (k)U° e (k, u) - V c (k)U° h (k, u) = l- V c (k) E^ln 

■ l 7T/C 

i=e,n 



u 2 - [(k 2 /2 mi ) + kv F , t 



(9) 



where vp, e /h is the Fermi velocity of electrons/holes at Fermi momentum in the conduction/valance 
band. 

In one-loop GW approximation with dynamically screened interaction, one has[8,33] 

z) = -- £ ^(^ *')^', z>) = -W -i^gL- Gi (A/, (10) 

^ fe',2' " fe',2' v ' / 

where 

and H e /h(k, z) is the self energy of electrons/holes defined in GW approximation. To avoid multi-pole 
structure in Gi(k,z) we approximate £j(fc, z) by momentum-dependent band gap renormalization 
Aj(fc). Using the approximation Aj(fc) = Ej(fc, — we find the following single pole electron- 
hole Green's function [8, 33] 

~r~T7T~ , (i = e,h) (12) 

2 — £i,k — Aj(/c) + //j 

The above formalism provides a suitable framework for our numerical calculation of band gap renor- 
malization. To proceed further, we should calculate the screened coulomb potential. Using equation 
(5), it can be written as follows 



V c (k) = ^ J dxdy J dx'dy'K \kR\\y(x,y)\ 2 \y(x',y')\ 2 , 
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or using rescaled coordinates Xi = ^ and x/ = ^ we find 



(13) 



where A is re-scaling factor equal to lCT 18 /m 2 and A;_R = k^J (x — rc ) 2 + (y — y ) 2 . Using ^(x,y) 
computed numerically in the previous section, we solve the integral of equation (13). The screened 
potential is calculated by ground state wave function as a function of k. Figure 15 shows the result 
of this calculation for different values of relative width of the confinement potential. In this figure, 
V c (k) is normalized by ^ and the k is normalized to kL. 

In the next step we calculate the band gap renormalization in both quasi static and GW approxima- 
tion. To do this end, we should calculate some quantities numerically. Using the re-scaled quantities 
k = kL and uo = Co x 10 16 , we write relation (9) as follows 



e(k, lo) — 1 




Co 2 x 10 7 



525[12fe - \ 



fel2 



Co 2 x 10 7 - 525[12|^ 



fcl2 



Co 2 x 10 fc 



174[43.1 



>k 2 
"1? 



fci2 
Li 



k]2 



(14) 
(15) 



x 10 6 - 174[43.8§ + L 

The single pole Greens function for electron defined as relation (12) transforms to the following form 

1 



G(k,z) = 



1.06 x 10- 26 iw 
and Greens function for hole becomes 

G(k,z) = 



6.17 x 10~ 21 g 



2.47 x 10- 19 - A e (k) 



(16) 



1.06 x 10- 26 i^ - 9.2 x 10- 20 || + 1.6 x 10~ 22 - A h (k) 



(17) 



To calculate bandgap renormalization in this configuration, we define the re-scaled j3 e and n e respec- 



tively as 



2m!W,2 



P 



574.5 



and i^P, where W y is width of quantum well in y direction in Nanometer. 



For holes, we also define the re-scaled f3h and //° as f3h — (3 e (^f) and n° h = ^(^f) respectively. In all 

h h 

computations in this paper, we have assumed that the ratio ^ is equal to 0.3 and m* ~ 0.067m e . 
Then we set Co = and using equation (7), we calculate band gap renormalization in both quasi static 
and GW approximations at temperature T = 0. The result of calculation for both quasi static and 
GW approximation is shown in figure 16. As this figure shows, the value of band gap renormalization 
in GW approximation is smaller than the quasi static plasmon-pole approximation. Typical values 
of band gap normalization are between 10-30 meV depending on the temperature and impurities 
in the system. Because of consideration of more quantum field theory details, GW approximation 
generally gives results which have better agreement with experimental results[10,ll]. 
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4 Interpretation of Numerical Results and conclusions 



Probability distribution of charge carriers in the absence of external electric and magnetic fields are 
shown in figure 2 where the asymmetric parameter a = ^ has been set equal to 0.8. This figure 
emphasizes the central role played by geometric shape of the confinement potential. In another words, 
in the absence of external electric and magnetic fields, carrier distribution obeys the symmetry of 
confinement potential. Variation of asymmetric parameter changes the profile of charge distribution 
in such a way that the case with a = 1 has maximum symmetry and any change of relative width 
leads to antisymmetric distribution. Figure 3 shows the variation of the ground state energy of 
carrier versus the inverse of the asymmetric parameter a. Variation of the relative width leads to the 
conclusion that smaller relative width leads to smaller ground state energy. Now, suppose that we 
turn on a uniform magnetic field. In the absence of electric field the distribution of charge carriers is 
given by figure 4. The role played by asymmetric parameter is a reduction of probability amplitude 
when the width of the well increases in x or y direction. Variation of ground state energy versus the 
intensity of magnetic field is depicted in figure 5 for two different values of asymmetric parameter. 
For a fixed well width in y direction, when the width of the well in x direction increases the value of 
the ground state energy will increase. 

Now we turn off magnetic field and apply a uniform electric field in the x direction. The profile 
of probability amplitude for charge carriers distribution is shown in figure 6 for a = 0.8. The 
probability amplitude has a Gaussian profile and is shifted toward the right hand side. This shift 
is a function of electric field intensity. The probability profile does not obey the geometric shape of 
the confinement potential due to preferred direction defined by the presence of electric field. The 
energy of the ground state versus the intensity of electric field is depicted in figure 7. As this figure 
shows, ground state energy decreases with increasing electric field intensity. This is not surprising 
since external electric field tends to decouple electrons and holes from each other. In the general 
case where both electric and magnetic fields are present, the asymmetric situation explained above 
will be enhanced in some respects. Figures 8 shows the space variation of the probability amplitude 
with a = 0.8. In the presence of constant electric and magnetic fields, pick of the graph describing 
charge carriers distribution will be shifted to the positive direction of the x axis. This feature causes 
the carriers concentration in such a way that these carriers distribution do not obeys the external 
confinement potential symmetry. More the intensity of the electric field results in more shift of 
the distribution pick to the right hand side of the x axis. The presence of constant magnetic field 
causes the anisotropy in the profile of the probability distribution of the carriers. Figure 9 shows 
the variation of the ground state energy versus the variation of the electric field when the magnetic 
field is supposed to be constant. More intensity of the electric field leads to the more reduction of 
the bound states energies. This resembles the linear Stark effect in elementary quantum mechanics. 
In the language of many body effects in dense plasma in a quasi one dimensional semiconductor, 
application of intense electric field results in the weaker excitonic states. Therefore, the presence of 
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external electric and magnetic fields will shift the location of maximum concentration of carriers and 
in this case there is an apparent asymmetry in profile of carriers distribution. This point can be used 
in fabrication of microelectronic devices based on quantum wires. 

We have proposed a numerical framework to calculate screened Coulomb potential and the values 
of band gap energy re-normalization in H-shaped quantum wire in two different approximations: 
quasi-static and dynamical random phase approximation in its leading order dynamical screening 
(GW approximation). We have evaluated the single particle self-energies for both electrons and holes 
in the dynamical plasmon pole approximation (PPA) and the leading order dynamically screened 
interaction or GW approximation to obtain the electron and hole renormalized Greens functions. This 
self-energy calculation gives us the band gap renormalization due to exchange-correlation effect. For 
comparison, we also calculated the band gap renormalization obtained by the quasi static calculation 
in both static random phase approximation and static plasmon pole approximation. Note that 
quasi static approximation works well in two and three dimensional systems but fails completely 
in one dimensional systems. This is because the electrons in one dimensional system suffer very 
strong inelastic scattering effects by virtue of restricted phase space. However, we have used this 
approximation only for comparison purposes. It is important to note that by re-scaling procedure 
which we have considered, we have fixed the geometry of the H-shaped quantum wire. Actually, one 
should consider the possibility for changing the geometry also. This has been down by change in 
parameter L but our calculations show that the main physical results do not change considerably. 

Figure 10 shows the screened Coulomb potential calculated based on random phase approxima- 
tion in its leading order and for different values of asymmetry parameter. Exchange-correlation many 
body effects mediate the bare Coulomb interaction. Based on different width of the confinement po- 
tential, screened Coulomb potential varies with geometrical characteristics of confinement potential. 
As the figure 10 shows, by increasing the asymmetry parameter, screened Coulomb potential will 
grows but its overall behaviors with respect to wave number will not change. Figure 10 shows the 
calculated band gap renormalization in quasi static and dynamical random phase approximations. 
To calculate band gap renormalization we first calculate the electron/hole single pole Green's func- 
tion and then using the formalism of both one-loop GW approximation and quasi static plasmon 
pole approximation we calculate the values of band gap renormalization in T = 0. Note that we 
have considered L dependence of band gap renormalization, however it can be translated to band 
gap renormalization versus carrier densities as well. Figure 11 also compares the results of band gap 
renormalization in GW and plasmon-pole approximations. As this figure shows, for a fixed value of 
L, GW approximation gives smaller values of band gap renormalization. Generally GW approxima- 
tion gives more reliable result in comparison with experimental data[5]. Note that we have focused 
on the variation of geometry by variation of L. As figure 1 shows there are other possibilities for 
changing the geometrical shape of external confinement potential. However the general behavior is 
the same as presented here. There are some restrictions on our calculations which can be summa- 
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rized as follows: the many-body treatment has the disadvantage that, for band gap renormalization 
it commonly ignores geometrical factors, such as the quantum confined Stark effect, whose relevance 
is structure specific. In other words, the numerical results for different confinement potential of H- 
shaped quantum wire may be geometry dependent in general. Furthermore, in the exciton problem, 
many-body theory treats screening within the linear approximation and, generally, influence of the 
bound electron on the free electrons is not fully included. In particular the orthogonality of the free 
electron states with the bound state, which increases its importance in lower dimensional systems, 
are typically not included. On the other hand the complete treatment of the problem should consider 
the effects of several subbands[8]. Note that in the rest of the calculation of band gap renormalization 
we have used the two band(one conduction band and one valance band) model and we have neglected 
the effects of higher subbands and degeneracies in valance bands. We also assumed that electrons 
and holes densities are constant in time. In summary, many body effects in quasi one dimensional 
systems are sensitive to the geometrical shape of external confinement effects. The value of band 
gap renormalization in GW approximation is smaller than quasi static plasmon-pole approximation. 
Typical values of this normalization are between 10-30 meV depending to temperature and impu- 
rities in the system. GW approximation generally gives results which have better agreement with 
experiment. A complete study of this problem requires the considerations of several subbands and 
varying electrons and holes densities. 



5 Summary and Conclusions 

Our numerical procedure to the bound states and many-body effects in H-shaped quantum wires has 
the following results: 

• The distribution of the probability of the charge carriers in the ground state of the H-shaped 
quantum wire in the absence of electric and magnetic fields has a symmetric shape obeying 
the geometric shape of the confinement potential. This distribution has its maximum at the 
center of each arm and decreases with distance from the center. The relative width of the 
confinement potential(the asymmetric ratio a = ^f- has considerable effect on the profile of 
this distribution. 

• In the presence of a constant magnetic field, the distribution of the charge carriers becomes 
oscillatory both in x and y directions. Increasing the strength of the magnetic field leads to the 
reduction of the ground state energy. The role played by the asymmetric parameter is given 
by the reduction of the probability amplitude when the width of the well increases in x or y 
directions. 

• The situation for the case of non-vanishing electric field (in the absence of magnetic field) 
resembles the Stark effect in a low dimensional system. The probability amplitude has an 
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oscillatory behavior with larger wavelength of the oscillations. In this case although the prob- 
ability amplitude has a Gaussian profile, it is shifted toward the right hand side. This shift is 
a function of the electric field intensity. The probability profile does not obey the symmetry 
of the geometric shape of the confinement potential due to preferred direction defined by the 
presence of the electric field. 

• In the presence of both electric and magnetic fields there are oscillations in probability distri- 
bution both in x and y directions, but in this case the probability distribution is not symmetric. 
In the presence of constant electric and magnetic fields, the pick of the probability amplitude 
will be shifted along the positive direction of the x axis. This causes the carriers to be concen- 
trated in such a way that they do not obey the external confinement potential symmetry. More 
intensity of the electric field results in more shift of the distribution pick to the right hand side 
of the x axis. The presence of constant magnetic field causes an anisotropy in the profile of the 
probability distribution of the carriers since it apparently breaks the local rotational symmetry 
in the center of each arm. 

• Screened Coulomb potential of the H-shaped external confinement is a sensitive function of the 
asymmetry parameter but its general behavior under variation of wave number is the same for 
other possible geometries of quasi one dimensional systems. The calculated band gap narrowing 
for H-shaped confinement potential in the absence of the electric and magnetic fields and within 
quasi-static and dynamical random phase approximation shows a typical gap narrowing of the 
order of few meV. This is supported from other studies of band gap narrowing in quasi one 
dimensional semiconductor systems[4,6]. The dynamical random phase approximation leads to 
more reliable result than quasi static approximation in comparison with experiments. 

In summary we can conclude that in the presence of electric and magnetic fields, H-shaped quan- 
tum wires bound states characteristics are sensitive functions of an asymmetric parameter a = |^ 
and the strength of the electric and magnetic fields. The case of non-vanishing electric and magnetic 
fields induces an intrinsic inhomogeneity in the quasi one dimensional system. Many-body effects 
due to plasma screening and resulting optical nonlinearities are also dependent to the asymmetric 
parameter of quasi-one dimensional confinement potential. Among these nonlinear optical effects, 
band gap renormalization of fundamental band age has been studied numerically in this paper. 

Appendix: Numerical Strategy 

We use the finite difference algorithm [34] to solve our partial differential equation(2) with boundary 
conditions imposed by H-shaped confinement potential. The most straightforward refinement method re- 
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places the differential equation with a finite difference equation. It replaces all derivatives with approximate 
expression as the following familiar form 



dy 
dx 



Vn+l ~ Vr, 
h 



dx 2 



Vn+l ~ tyn + Vn-1 
h 2 



(18) 



where the mesh chosen to be equally spaced and given by h = ^+1° with xt = XQ+kh and k = 1, 2, n+1. 
We first re-scale the Schrodinger equation and the confinement potential. We define the re-scaled value of x 
and y as x = and y = where W x and W y are the width of the well in x and y direction respectively. 
Now the re-scaled Schrodinger equation can be written as follows 



i d 2 



1 d 2 2ieBW x _d\ eBW x 

\w x dx 2 + w y dy 2 hw y x dy) (x,y)+ { h j 

2meEW x _ 2m 

to x^{x,y) = -2-A^(a;,y). 

n n 

Using the finite difference algorithm, this equation can be written as follows 

( - 2* hj + o jW " + 

~T~ CI 



x 2 ^(x, y) 



A, 



(19) 



+ 



UeaBW* _ * <J+1 - / eBW 2 



2A, 



+ 



n 



2meEWl 



n 2 



x x&ij = 



2mWl 



(20) 



where a 



Wy 



is the relative width of the H-shaped quantum wire. We do discritize the x and y axes to 



discrete space dx = 0.1 and dy = 0.1, therefore the equation (20) can be written as a matrix equation 

H<p = E<f>, (21) 
where Tl, and <j) are the Hamiltonian matrix and the state wave functions array which are defined as follows 



/ 



\ 



4>(xi,yj) 

(j)(x i+1 ,yj) 

V : ) 

4>(xi,y j+ i) 
<f>(xi+i,y j+1 ) 



(22) 
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V 



K 



. \ 





-) + (^) 2 



3 



2 (^T + - 



V <> 



V o 



(23) 



The error of computation of 4>(x,y) is of the order of 0(dx 2 ). We diagonalize the Hamiltonian 
and calculate the bound states energies and wave functions numerically. We use MATLAB package 
uses techniques that are more efficient than Jacobi rotations and that can be applied to asymmetric 
complex matrices as well as to the more common real symmetric situations. 
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Figure 1: The Geometry of the H-Shaped Quantum Wire 
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Figure 3: The ground state energy versus the inverse of the asymmetry parameter at zero external fields 
strength. 
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Figure 4: The probability amplitude of charge carriers distribution for ground state of H-shaped quantum 
wire in the presence of an external magnetic field(and zero electric field). The asymmetry parameter has 
been chosen to be a = Wr = 0.8. 



18 



9.5 1 

9-0*, 
8.5 - 

8 - 
7.5 - 

7 - 
6.5 - 

6 - 
5.5 - 

5 - 



■ -O- i a=1 .2 
>-•-•' a=1 



"0^ \ 



V, 



15 



20 



25 



30 



3.06 
3.04 
8.02 

I 

8 

7.98 
■7.96 
■ 7.94 



^7.92 
35 



Figure 5: The ground state energy versus the intensity of the external magnetic field for different values of 
asymmetry parameter. 
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Figure 6: The probability amplitude of charge carriers distribution for ground state of H-shaped quantum 
wire in the presence of an external electric field(and zero magnetic field). The asymmetry parameter has 
been chosen to be a = 0.8. 
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Figure 7: The ground state energy versus the intensity of the external electric field for a = 0.8 and B = 0. 
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Figure 8: The probability amplitude of charge carriers distribution for ground state of H-shaped quantum 
wire in the presence of both electric and magnetic fields. The asymmetry parameter has been chosen to be 
a = 0.8. 
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Figure 9: The ground state energy versus the intensity of the external electric field for a = 0.8 and B = 20 
(a constant magnetic field). 
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Figure 10: The calculated screened Coulomb potential versus the wave number for different asymmetry 
parameter. The screened potential is normalized by ^— and the k is normalized by L = W y . 
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Figure 11: Calculated band gap renormalization in H-Shaped quantum wire within dynamical random 
phase(GW) approximation (upper curve) and quasi-static plasmon-pole approximation (lower curve). The 
screened potential is normalized by — and the k is normalized by L = W y . 
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